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Abstract 

We use an infinite-range Maier-Saupe model, with two sets of lo- 
cal quadrupolar variables and restricted orientations, to investigate 
the global phase diagram of a coupled system of two nematic subsys- 
tems. The free energy and the equations of state are exactly calculated 
by standard techniques of statistical mechanics. The nematic- isotropic 
transition temperature of system A increases with both the interaction 
energy among mesogens of system B, and the two-subsystem coupling 
J. This enhancement of the nematic phase is manifested in a global 
phase diagram in terms of the interaction parameters and the tem- 
perature T. We make some comments on the connections of these 
results with experimental findings for a system of diluted ferroelectric 
nanoparticles embedded in a nematic liquid-crystalline environment. 

1 Introduction 

A dilute suspension of ferroelectric nanoparticles in a liquid-crystalline host 
has been shown to display an enhancement of the nematic order, with an 
increase of the isotropic-nematic transition temperature, which is a behavior 
of interest from the point of view of technological applications (TJ, [2j [3] . This 
effect has been explained by the introduction of a coupling between the usual 
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nematic order parameter of the liquid crystals and a set of extra degrees 
of freedom associated with a coarse-graining average of the electric dipole 
field produced by the ferroelectric nanoparticles jU |5]. This work was the 
motivation to investigate the global phase diagram of a basic Maier-Saupe 
(MS) model on a lattice with two coupled sets of quadrupolar degrees of 
freedom. 

Although the nematic-isotropic transition is perhaps the most explored 
transition in liquid crystalline systems, there are still a number of questions 
and open problems, which can be formulated in terms of simple statistical 
lattice models. An interesting question is the onset of a biaxial nematic 
phase [5], which we have recently investigated in the context of a MS model 
for a mixture of discs and cylinders [7J [5]. Now we analyze the global phase 
diagram of a similar type of statistical model, with the inclusion of two sets 
of quadrupolar degrees of freedom, which leads to a connection with the 
work by Lopatina and Selinger [U [5]. The nematic MS model is the liquid- 
crystalline analog of the Curie- Weiss model of ferromagnetism [9], [10l ITT] . 
In this approach, the standard nearest-neighbor interactions between lattice 
sites are adequately replaced by scaled interactions of infinite range. The 
statistical-mechanics problem is exactly solvable, and leads to a very con- 
venient framework to perform calculations at the mean-field level. The MS 
model can be further simplified if we suppose that the local mesogen orienta- 
tions are restricted to a discrete set of states, according to an early proposal 
by Zwanzig |12j . Some model calculations with continuous orientations for 
uniaxial [13] . and biaxial [8] nematic systems give support to the idea that 
this simplification does not lead to qualitatively different results. In recent 
publications, we have used extensions of this Maier-Saupe-Zwanzig (MSZ) 
lattice model to investigate the existence of biaxial nematic phases [HI [71 [J3] 
and the thermodynamic properties of nematic elastomers [TS] . 

In Section II we define the MSZ model with two sets of coupled degrees of 
freedom, and use standard tools of statistical mechanics to write a thermo- 
dynamic free energy. This approach has a number of advantages. In contrast 
with the standard Landau phenomenological approach, the microscopic in- 
teractions are explicitly stated and the calculations are not restricted to the 
neighborhood of the transitions. Also, the model is simple enough to allow 
exact bona-fide calculations of the free energy and equations of state. In 
Section III we study the global phase diagram. Contact with a dilute system 
of ferroelectric nanoparticles embedded in a nematic host is made in Section 
IV. We then conclude with a summary of the main results, which do support 



2 



the experimental enhancement of the nematic order. 



2 Coupled MSZ model 

The energy of a model with two coupled subsystems A and B can be written 

as 

n = E A + E B + E AB , (1) 
where E\ and are given by standard Maier-Saupe soft quadrupolar forms, 

^a = -| E E ^r, (2) 

l<i<j<N fi,v=x,y,z 

Eb = E E Brer, (3) 

l<i<j<N fi,u=x,y,z 

with the quadrupole components 

at = \ (3 af < - <n , er = ^ (3 - <n , (4) 

where {a,j} and {6j} are the respective sets of unit vectors associated with 
mesogenic units of types A and B, e A and 5b are energy parameters, 5^ v is a 
Kronecker symbol, and the coupling term is given by 

£ab = ~ e e A^r. (5) 

l<i,j<N n,v=x,y,z 

where J is the coupling parameter. Note that the energy global forms for 
E&, Eb, and E AB are consistent with the mean- field level of calculations 
considered in this paper. 

In the MSZ model, calculations of the canonical partition function involve 
sums over the local orientations of the mesogenic units, which are restricted 
to the six directions along the three Cartesian axes, 

a t , b % e {(±1, 0, 0), (0, ±1, 0), (0, 0, ±1)}. (6) 

The partition function can be written as 



z = EE ex p^E 

{at} {bi} I n,v 



N \ / N 



+ 2j eah E^ 



j=i j \i=i 
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where (3 = l/(kT) is the inverse temperature, k is the Boltzmann constant, 
and we have discarded irrelevant terms in the thermodynamic limit. In order 
to linearize the quadratic forms, we consider the set of identities 



nw) n 



fJ,,U 



fJ,,U 
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2~Ttl 
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E4" 



NQ>? 
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N 
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1, (9) 



where we have used the complex representation of the Dirac 5-function. We 
then write 



exp 



[dz] 
/3N 



E 



fl,U 



[dQ] 
10) 



where [dz] = [] dz^dz^ 1 , [dQ] = Y\ v dQ^dQ 1 ^, and we have discarded 
contributions to the free energy of 0(lniV). Since the interaction energies, 
given by Eqs. 02H5]), are of infinite range, the sum over states has been 
reduced to a simple problem, with decoupled sites, and the final results are 
of mean-field nature. The sum over the A variables leads to 



E ( E E # V) = E (e *aaA 



N 



exp 



ln2-l^,r + ln [£, 



N 



(11) 



which also holds with the exchange A-H-B. 

For large N, the integrals over the Q variables can be obtained by the 
Laplace method. To leading order, we have 



Iq = exp 



N 



2/3 (e A e B - J 2 ) 



(12) 
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The stationary point conditions, 



z 



P ( £ aQ*a + JQ 



B I > 



% = + j Qa 

which can also be written as 



A " /J {£a£b _ J2) > p {£a£b _ J2) 



(13) 



(14) 



may be used to switch to the former Q variables. Eqs. (fTTj) and f lT2|) can be 
used to write the partition function 



27V In 2 



[dz] exp I -N 



2/3 (e A e B -Ji) 



+ \ £ (4* + *T) - In ( £ e 3 ^ 2 ) - In ( £ 



,3^/2 



Note that, on the grounds of mathematical simplicity, we decided to begin 
the integrations by the Q rather than the z variables. 

The complex integrals over the z variables may be done by the method 
of steepest descents. For \i ^ u, the intergral over z^ v gives a contribution of 
O(lniV), which will then be discarded. The saddle-point conditions for the 
diagonal variables lead to the self-consistent equations 



or 



Me A Q^+JQ>^)/2 



3f3(e A QT+JQB)/ 2 



(16) 



(15) 



Q 



B 



J2 e 3/3(e B Q^+JQr)/ 2 



(17) 



where we have switched to the Q variables using Eqs. (fTBl and ( !T4l) . Note 
that we recover the simple MSZ model in the limit J — y 0, and that the Q 
tensors are traceless. The free energy is finally written as 

/ = -2fcTln2 + \ ( £ aQT 2 + e B Q B v2 + 2JQ K V Q™) - kT 



In 



£exp (^p( £A QT + JQ^ 



+ In 



J2exp(-f3(e B Q£ + JQ A n 
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3 Global phase diagram 



The equations of state and the free energy are further simplified if we consider 
the standard parametric form of the order parameter tensors, 



(19) 



5 5 a 
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S A 


5-Sb 











-§<Sb 











Sb 



(20) 



so that 



5 A 
Sb 



1 + 2 e -9(SA+JS B )/*T 

3 



1 + 2 e -9(e B SB+JS A )/4T 



- 1 
- 1 



(21) 



and 
/ = 



-2Tln2 + - [S A (S A - 2) + BbSb(S b - 2) + 2J(S A S B - (3 A + S B ))} 

(22) 



-Tin [1 + 2 e _9(5A+J5B)/4T ] [l + 2 e" 9(£B5B+J5A)/4T ] . 



Note that we have written the free energy / and the energy parameters Eb 
and J, as well as the temperature T, in terms of the energy parameter e A . 

Assuming eb < s A , we can write an expansion of the free energy / in 
terms of S A , up to 0(S A 2 ), 

f = -2 Tin 6 + j_ ( 3g B~3J 2 -4T) _ 2 

16 (3e B -4T) 2 T L 1 A ;j a , i ; 

so that the coefficient of S"a 2 is proportional to T — T*, where T* is given by 

T* = (l + e B + v / (1+^b) 2 + 4(J 2 -£ B )) • (24) 

If we suppose that the coefficients of higher powers of S A in f )23|) are weakly 
dependent on temperature, T* is associated with the spinodal temperature 
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of a Landau-de Gennes theory of nematics, and establishes the threshold of 
stability of the isotropic phase. For e b — > 0, the spinodal temperature is 
given by T* = 1 + VI + J 2 , which increases with the coupling parameter 
J, and is thus an indication of the enhancement of the nematic phase. For 
sb — > 1, T* increases approximately linearly with J. The global behavior is 
shown in the diagram of Fig. [U in terms of Eb, J, and T*. In particular, 
the coefficient of S\ 2 becomes proportional to (T — 3/4) in the limit of zero 
coupling between the subsystems, J — > 0, in agreement with previous results 
for the simple Maier-Saupe-Zwanzig model [7J. 




Figure 1: Stability threshold of the isotropic phase in terms of the spinodal 
temperature T* and the energy parameters Eb and J. The isotropic phase is 
unstable below the surface. 

In Fig [21 we show the global phase diagram of this model in terms of the 
temperature T and the energy parameters Eb and J, obtained from numerical 
solutions of the exact equations of state ( 12TI) and free energy ( |22i) . The transi- 
tion from the nematic (low-temperature) to the isotropic (high-temperature) 
phase is of first order, with a jump in both nematic order parameters S\ and 
Sb- The nematic phase is enhanced by the increase of the energy parameters 
Eb and J. As it should be anticipated, this surface is very similar to the 
threshold of stability of the isotropic phase. 
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Figure 2: Global phase diagram in terms of temperature T and energy pa- 
rameters £b and J. The transition from the nematic (below the surface) to 
the isotropic phase (above the surface) is discontinuous, with a jump in the 
nematic order parameters Sa and Sb- 

4 Suspensions of ferroelectric nanoparticles 
in nematic systems 

A coupled system of two types of nematic subsystems has been considered 
by Lopatina and Selinger to represent a dilute suspension of ferroelectric 
nanoparticles in a nematic host [U |5] . The mechanism behind this mapping 
is the effect on the nematic mesogenic units of the dipole aligning electric 
fields produced by the ferroelectric nanoparticles. For a range of param- 
eters, assuming uniform nematic order, one may integrate out the position 
variables, and thus eliminate the complicated spatial dependence of the inter- 
actions between mesogens and nanoparticles. It is then possible to associate 
a nematic-like order parameter with the distribution of orientations of the 
dipole moments of the nanoparticles. The energy of interaction turns out 
to be proportional to Slc ' ^npj where Slc and Snp represent the nematic 
order parameters of the liquid crystal and nanoparticle systems respectively 
[U [5] , and which has essentially the same form as the coupling energy of Eq. 
dSJ). Though the interaction among the nanoparticles has dipole symmetry, 
one expects the limit eb — > will approximately represent a situation in the 
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dilute regime. 




Figure 3: (a) Nematic order solutions for S\ (full curves) and Sb (dashed 
curve) as a function of temperature for e B = 0, and some values of J. There 
may be up to three solutions: stable nematic, unstable nematic, and isotropic 
(S = 0). The transition from the nematic (upper solution) to the isotropic 
phase (5* = 0) takes place where the two free energies become equal, (b) Free 
energy as a function of temperature for the stable and unstable nematic (full 
curves) and the isotropic solution (dashed curve), with e B = and J = 0.5. 
The inset shows a magnification of the region where the stable nematic and 
the isotropic solutions have the same free energy. 

In Fig [3h we show the nematic solutions for S\ (full curves) and 5b 
(dashed curves) as a function of temperature, for e B = and some values of 
the coupling parameter J. In general, there may be up to three solutions, 
two nematic and one isotropic. The lower curve is unstable, and the nematic- 
isotropic transition takes place at the temperature at which the stable ne- 
matic (upper curve) and isotropic (S = 0) free energies become equal. The 
free energies for the three solutions are shown in Fig [3b, for J = 0.5. In agree- 
ment with some previous experimental [2] and theoretical studies [H [3, [3] , 
the nematic order is enhanced by the coupling with an ordered subsystem 
B, which may be interpreted as representing the set of local orientations of 
ferroelectric nanoparticles. Note that our model always predicts an increase 
of Tni with the addition of ferroelectric nanoparticles, even though a scenario 
where the nematic-isotropic transition temperature decreases has also been 
observed [17]. This scenario is probably suppressed by the assumption of 
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uniform nematic order around the nanoparticles, which our model inherits 
from Lopatina and Selinger's derivation of the nematic-nanoparticle coupling 
Hi- 



5 Summary 

We have used a Maier-Saupe-Zwanzig lattice model to study a coupled sys- 
tem of two types of nematic subsystems. The model is simple enough to 
allow an exact calculation of the free energy. In contrast to the usual phe- 
nomenological approaches, it explicitly contains the microscopic energy pa- 
rameters involved in the nematic-isotropic transition. We show that the 
nematic-isotropic transition temperature of subsystem A increases with both 
the interaction energy among mesogens of system B and the two-subsystem 
coupling J. We draw the global phase diagram in terms of the temperature T 
and the energy parameters. For non- interacting mesogens in system B, this 
model may be used to describe the experimentally observed enhancement of 
the nematic ordering produced by ferroelectric nanoparticles suspended in a 
liquid-crystalline host. 
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